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Abstract. The convergence of solutions of the incompressible 
Navier-Stokes equations set in a domain with boundary to solu- 
tions of the Euler equations in the large Reynolds number limit is 
a challenging open problem both in 2 and 3 space dimensions. In 
particular it is distinct from the question of existence in the large 
of a smooth solution of the initial-boundary value problem for the 
Euler equations. The present paper proposes three results in that 
direction. First, if the solutions of the Navier-Stokes equations sat- 
isfy a slip boundary condition with vanishing slip coefficient in the 
large Reynolds number limit, we show by an energy method that 
they converge to the classical solution of the Euler equations on its 
time interval of existence. Next we show that the incompressible 
Navier-Stokes limit of the Boltzmann equation with Maxwell's ac- 
commodation condition at the boundary is governed by the Navier- 
Stokes equations with slip boundary condition, and we express the 
slip coefficient at the fluid level in terms of the accommodation 
parameter at the kinetic level. This second result is formal, in the 
style of [Bardos-Golse-Levermore, J. Stat. Phys. 63 (1991), 323- 
344]. Finally, we establish the incompressible Euler limit of the 
Boltzmann equation set in a domain with boundary with Maxwell's 
accommodation condition assuming that the accommodation pa- 
rameter is small enough in terms of the Knudsen number. Our 
proof uses the relative entropy method following closely [L. Saint- 
Raymond, Arch. Ration. Mech. Anal. 166 (2003), 47-80] in the 
case of the 3-torus, except for the boundary terms, which require 
special treatment. 
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I. Introduction 

In a program initiated more than 20 years ago with Dave Lever- 
more [3, 4, 5], we outlined a strategy for deriving incompressible fluid 
dynamic equations from the theory of renormalized solutions of the 
Boltzmann equation invented by R. DiPerna and P.-L. Lions [12]. 

At the time of this writing, complete derivations of the Stokes [24, 
25], Stokes- Fourier [14] and Navier- Stokes- Fourier [15, 16, 20] have been 
obtained following that program, in the greatest possible generality 
allowed by the current existence theories for both the fluid dynamic 
and the Boltzmann equations: see [36] for a survey of these issues. 

The case of the incompressible Euler equations in space dimension 3 
stands out, in the first place because there does not exist a satisfactory 
theory of global weak solutions of these equations analogous to Leray's 
theory of weak solutions of the Navier-Stokes equations [21] in space di- 
mension 3. Even if there was a global existence theory of weak solutions 
of the incompressible Euler equations in the energy space L^{L 2 X ) in 
dimension 3, such solutions would not satisfy the weak-strong unique- 
ness property observed by Leray in the case of the Navier-Stokes equa- 
tions. (Indeed there exist nontrivial compactly supported solutions of 
the incompressible Euler equations in energy space: see [32, 33, 11].) 
In [23], P.-L. Lions proposed a notion of dissipative solution of the 
incompressible Euler equations — in the same spirit of his definition 
of the notion of viscosity solutions of Hamilton- Jacobi equations, but 
using the conservation of energy instead of the maximum principle as 
in the Hamilton- Jacobi case. The weak-strong uniqueness property is 
verified by dissipative solutions of the incompressible Euler equations 
(essentially by definition): if there exists a classical (C 1 ) solution of 
the incompressible Euler equations, all dissipative solutions with the 
same initial data must coincide with this classical solution on its maxi- 
mal time interval of existence. Unfortunately, dissipative solutions are 
not known to satisfy the incompressible Euler equations in the sense of 
distributions. 

Using the relative entropy method pioneered in [38] and adapted 
to the case of the Boltzmann equation in [10, 25], L. Saint-Raymond 
[28, 29] succeeded in deriving dissipative solutions of the Euler equa- 
tions in arbitrary space dimension (or classical solutions whenever they 
exist) from weak solutions of the BGK model [28] or from renormalized 
solutions of the Boltzmann equation [29]. 

However, all the derivations of fluid dynamic equations from the 
Boltzmann equation referred to above are carried out in either the 
Euclidian space of the flat torus so as to avoid difficulties 
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that may result from boundary conditions. The theory of renormalized 
solutions of the Boltzmann equation in the presence of accommodation 
boundary conditions was obtained only very recently, by S. Mischler 
[27]; subsequently, N. Masmoudi and L. Saint-Raymond established the 
Stokes- Fourier limit of such solutions [26]. 

In the present paper, we derive dissipative solutions (or classical 
solutions whenever they exist) of the incompressible Euler equations 
from renormalized solutions of the Boltzmann equation in some spatial 
domain satisfying Maxwell's accommodation boundary condition. In 
particular, we identify a sufficient scaling condition on the accommo- 
dation parameter under which the hydrodynamic limit of the family of 
solutions of the Boltzmann equation is governed by the incompressible 
Euler equation with its classical boundary condition — i.e. assuming 
that the velocity field is tangent at the boundary. 

The outline of the paper is as follows. Section 2 gives a sufficient 
condition on the slip coefficient at the boundary under which the in- 
compressible Euler equations are obtained as the inviscid limit of the 
incompressible Navier-Stokes equations with slip-boundary condition. 
The main result in this section is Theorem 2.2, based on an energy 
method. Section 3 provides a formal derivation of the incompress- 
ible Navier-Stokes equations with slip boundary condition from the 
Boltzmann equation with Maxwell's accommodation condition at the 
boundary of the spatial domain: see Theorem 3.1 for a precise state- 
ment of this result. Based on the intuition provided by sections 2 and 
3, we identify a scaling limit of the Boltzmann equation with Maxwell 
accommodation boundary condition leading to the incompressible Eu- 
ler equations: see Theorem 4.2, whose proof occupies most of section 
4. 

It is a our great pleasure to offer this modest contribution to our 
friend Dave Levermore, in recognition of his outstanding influence on 
the analysis of nonlinear partial differential equations in the past 30 
years, especially on the problem of hydrodynamic limits of the Boltz- 
mann equation, directly inspired from Hilbert's 6th problem on the 
axiomatization of physics. 



2. Inviscid Limit of the Navier-Stokes Equations 
with Slip Boundary Conditions 

A a warm-up, we begin with a simple observation bearing on the 
inviscid limit of the incompressible Navier-Stokes equations set in some 
smooth domain with slip boundary condition. In particular, we identify 
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a sufficient scaling condition on the slip coefficient in order to obtain 
the incompressible Euler equations in the inviscid limit. 

Let Q designate an open set in H N with C 1 boundary dfl, assuming 
that N = 2 or 3; henceforth the outward unit normal vector at the 
point x of dfl is denoted by n x . Consider the initial-boundary value 
problem with unknown u v = u u (t, x), set for x G Q and t > 0: 

( div x u v = , 
d t u v + dw x (u u ® u v ) + V x p v = v& x u v , 



(1) 



u u -n\ m = 0, 

v{H{u v ) ■ n) T + AwJ = , 



Uu\ t=Q = U m 



where v > is the kinematic viscosity, A > the slip coefficient, 
E(w) := V ' x u + (V x u) T , while 

v(x) T := (I — n(x)® 2 ) ■ v(x) . 

Henceforth, we denote 

U(n) := {v G L 2 {VL; R N ) | divv = and v ■ n\ m = 0} . 

For each v > and -u m G "H(fi), there exists a weak solution u v of (1) 
in L°°(R + ; fl L 2 (R + ; meaning that, for each test vector 

field [/ G C(R + ;'H(Q)) n C C °°(R+ x H), one has 
(2) 

/ : Y,{U)dxdt + X f I u v ■ UdS(x)dt 

Jn Jo Jan 



roo r f 

I I iu v ■ d t U + u u ®u u : V x U)dxdt + / - 
Jo Jn Jn 



u in {x) ■ U(0,x)dx, 



and satisfying in addition u u G C(R + ;-u; — L 2 {Vt)), together with the 
Leray-type energy dissipation inequality: 

/ 7;\u u (t, x)\ 2 dx + v / / \E(u v )(t, x)\ 2 dxdt 
^ Jn Jo Jn 

+ X I [ \u u (t,x)\ 2 dS(x)dt< [ \\u m {x)\ 2 dx 
Jo Jdn Jn 

for each t > 0. Such a weak solution of (1) will henceforth be referred 
to as a "Leray solution of (1)". The classical theory of Leray solutions 
that is well known in the case where the velocity field satisfies the 
Dirichlet boundary condition on dQ can be adapted to the case of the 
slip-boundary condition: see [34, 7], and Theorem 2 in [18]. 
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Since u v G C(R + ;w — L 2 (Q)), by an elementary density argument, 
one can choose a sequence of test vector fields U n of the special form 
U n (t,x) = Xn(T-t)w(t,x) where w G C(R+; H(n)) nC c 1 (R + x H) and 



Xn(z) = / Xn(s)d.S 
J -co 

where x' n * s a regularizing sequence on R, so that the weak formulation 
(2) of the Navier-Stokes equations becomes 

v \ / |E(u„) : Yu{w)dxdt + X / / u v ■ wdS(x)dt 
Jo Jn Jo JdQ. 

u m {x) ■ w(0, x)dx — / u u (T, x)w(T, x)dx 
n Jn 

r T 



/ / (u v ■ d t w + u v ® u v : V x w)dxdt , 
Jo Jn 



for each T > 0. Furthermore, denoting 

E(w) := d t w + w ■ V x w , 



one has 

rT 



' / u v ■ d t wdxdt — / [u v ■ E(w) — u u ® w ■ V x w)dxdt 
o Jn Jo Jn 



while 

rT 



\ \\w{t,x)\ 2 dx — I ||w(0, x)\ 2 dx + / / w ® w :V x wdxdt 
Jo Jo Jo Jn 



(4) ^ 



Therefore 



o Jn 

/ w ■ E(w)dxdt . 
o Jn 



v 

rT 



I I \Ti{u v ) :Ti{w)dxdt + X / / u u -wdS(x)dt 
'o Jn Jo J on 

(5) = / / [u v ■ E{w) + u v <8> — iu) : V x w — w ® u p : V x w)dxdt 
Jo Jn 

+ / u m (x) ■ U(0,x)dx- / u v (T, x)w(T, x)dx , 
Jn Jn 



since 



Jn 
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because u v • n\ BQ — 0. Combining (5) and (4), we find that 



(6) 



v 



JCl 



I / |E(u„) : Ti{w)dxdt + A / / u u ■ wdS(x)dt 
o Jn Jo Jon 

{{u v — w) ■ E{w) + [u v — w) <S> (u u — w) : V ' x w)dxdt 

+ / u in (x) ■ w(0,x)dx- I ±\w{0,x)\ 2 dx 
Jn Jn 

I u u (T, x)w{T 1 x)dx + / \\w(T, x)\ 2 dx 
n Jn 



Finally, combining (3) and (6), we conclude that any Leray solution u v 
of (1) satisfies the inequality 



/ \\u u — w\ 2 {t,x)dx + I I (u v — w) ® (u v — w) :V x wdxds 
Jn Jo Jn 



(7) 



+v / / ||E(u I/ )(s, x)\ 2 dxds + A / / \u u (s, x)\ 2 dS(x)ds 
Jo Jn Jo Jon 

< [ \\u ir \x) - w(0,x)\ 2 dx+ [ [ E(w) ■ {u v - w)dxdt 
Jn Jo Jn 

+u / / \Y,{u v ) : T*(w)dxds + A / / u v • wdS(x)ds 
Jo Jn Jo Jon 



for each w e C(K+; ft(fi)) n C C 1 (R+ x Q). 

At this point we recall the definition of dissipative solutions of the 
incompressible Euler equations set in a domain Q with smooth bound- 
ary: 



(8) 



diva; u — , 

d t u + &iv x (u ®u) + V x p = , 



u ■ n\ 



an 



0, 



u 



t=0 



u 



Definition 2.1 (P.-L. Lions [23], p. 154, C. Bardos, E. Titi [6], p. 
16). Given u m G H{p), a dissipative solution of (8) is an element 
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u G C(R + ; w — "H(fi)) satisfying u 



t=o 



u m and the inequality 




u — w\ 2 (t, x)dx 




\Jo J Jn 



+ / exp I / 2 1| cr(w) \\L°°(n)( r )dT J / E{w) ■ (u — w)(s, x)dxds 




for each w G C(R+; H(£l)) n C\R + x Q), where 



(9) 



a(w) (t,x) : = sup(-S(w)(t,x) : £ ® £) . 
I«l=i 



We recall that, if the Euler equations (8) have a classical solution 
v G C^QO.T*) x H) satisfying 

<t(u) _ GL^O,!];! 00 ^)) and p G ^([0, T]; tf 1 ^)) for each T < T* , 

then all dissipative solutions of (8) must coincide with v on [0, T*) x f2 
a.e., since one can use w = v as the test vector field, so that 



because (u — v)(s, •) G "H(f2) for each s G [0,T). 

Theorem 2.2. Let u m G "H(fi), and assume that the slip coefficient 
A = \{u) in (1) scales with the kinematic viscosity u so that 



Then any family (u u ) of Leray solutions of (1) is relatively compact in 
the weak-* topology o/L°°(R+; and in C(R+; w - for the 

topology of uniform convergence on bounded time intervals, and each 
limit point of (u u ) as v — > is a dissipative solution of (8). 

Proof. We deduce from (7) with w = 0, or equivalently from the Leray 
energy inequality that 





(10) 



\(u) -> as v — >■ . 



(11) 



\JvYi(u v ) is bounded in L 2 (R + ; L 2 (f2)) , and 



^/\(u)u u \ dn is bounded in L 2 (R + ; L 2 (dQ)) . 



8 C. BARDOS, F. GOLSE, AND L. PAILLARD 

By Gronwall's inequality 
(12) 

/ \\u u — w\ 2 (t, x)dx 
Jn 

<exp^y 2\\a(wy\\ L oc {n) (s)ds S j J \\u in (x) - w(0, x)\ 2 dx 

+ / exp I / 2\\a(wy\\ Lx ^- ) (r)dr) / E(w) ■ (u v - w)(s,x)dxds 
Jo \Js J Jn 

+ J exp ^ 2\\a(wy\\ L o ^ n) (T)dr^ Q u (s)ds 

where, by the Cauchy-Schwarz inequality 

Q v {s) = u\\T,(u u )\\ L 2 {n) (s)\\T,(w)\\ L 2 {n) (s) 
+ Hv)\\ u v\\L2(dn)(s)\\w\\ L 2 {n) (s) . 

In view of (11), one has 

HaiUw]) = O(v^) + 0{y/My)) -+ 

as v — > 0, and we conclude by passing to the limit in (12) following the 
same argument as in [23]. □ 

Several remarks are in order after this result. 

In some references, the slip boundary condition is written 

(13) u(^) +\u\ dn = 0, 



dn / 



instead of 

(14) u(E(u v )-n) T + Xu\ Xi = 0. 

Likewise, the boundary condition 



curl u v x n\ m = 

is also considered by some authors — and referred to as the Navier slip 
condition — in the context of the inviscid limit of the Navier-Stokes 
equations: see for instance [22, 2, 7, 9]. 

If dQ is a straight line, or a plane, or a hyperplane in space dimension 
N > 3, the normal vector field n is constant, so that 

(£(«,) -n)J= -f +V>-nL Q ) = 



an \ dn I l5n \dn J \ dn I l5n 

"' T 
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and 



curlw,, x n\ m = {V x u v - (V ' x uf) ■ n\ m 

du v \ , t . i I 'du 



g n y I an v \an> \ dn I 9n 

(with V T denoting the tangential component of the V operator), since 
the velocity field u v is tangential on dVt. Therefore, in the case of a flat 
boundary, all these boundary conditions are equivalent. 

If dVL is a smooth curve, or a surface or a hypersurface in space 
dimension N > 3, then 

(T,(u u ) ■ n) T \ = ( ^ J -V'ra-u 



9n \ dn ) T \dfi 



curl -u^ x n\ „„ = I ^— - I + V T n • u 1 



while 

l9n \ dn ) ' r lan 

so that all these boundary conditions differ by a 0-order operator given 
by the Weingarten endomorphism of the boundary dQ. 

Here, we have chosen the second boundary condition above, as it is 
the more natural one when looking at the Navier-Stokes equation as a 
fluid dynamic limit of the kinetic theory of gases. 

However, the same argument as in the proof of Theorem 2.2 can be 
extended to treat the case of a slip coefficient A which is not nonnega- 
tive, provided that 

A(z/) + = max(A(z/),0) and A_(i/) = max(-A(z/), 0) = 0{u) 

as v — y 0. Indeed, the contribution of A(z/)_ in the estimate (7) can 
be absorbed in the viscous dissipation term by means of the following 
classical inequality: for each a > 0, there exists C a > such that, for 
each v G H^Q), 



dS(x)<a [ \V x v(x)\ 2 dx + — [ \v(x)\ 2 dx. 
Jn « Jq 



>dfl 

With this observation, the term 

(£(u„) • n 

can be replaced indifferently with either 

' du v 



an 



dn 



an or cur\u u xn\ m 



in the slip boundary condition. 

More precise variants of Theorem 2.2 have been established by var- 
ious authors, see for instance [37, 8]. The result given here holds for 
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a very general class of nonnegative slip boundary coefficients A and is 
based upon the simplest imaginable energy estimate. The condition 
A > in (1) is somehow natural when this initial-boundary value prob- 
lem is considered as some scaling limit of the Boltzmann equation of 
the kinetic theory of gases. 

Another question is whether the condition X(u) — y as v — » is 
optimal. Considers instead the Dirichlet boundary condition for u u , 
i.e. 

UL, \an = • 

Formally, this boundary condition corresponds with any one of the slip 
boundary conditions above with 

lim > 0. 

In that case, it well known that the Euler equations (8) may fail to 
describe the inviscid limit of the Navier-Stokes equations, even in the 
simpler 2 dimensional case. Because the Dirichlet boundary condi- 
tion overdetermines the velocity field in the inviscid limit, the Euler 
equations (8) are expected to govern the inviscid limit of the Navier- 
Stokes equations only if the effect of viscosity remains confined on a 
thin layer near the boundary. But it may happen — and does happen 
under certain circumstances — that the viscous layer detaches from 
the boundary, as for instance in the case of the so-called von Karman 
vortex streets in the case of a Navier-Stokes flow past a cylinder, even 
at moderate Reynolds numbers. While this situation seems beyond 
the grasp of current mathematical analysis, there exists a least a very 
interesting criterion due to T. Kato [19], formulated in terms of the 
viscous energy dissipation only, identifying situations where the invis- 
cid limit of the incompressible Navier-Stokes equations with Dirichlet 
boundary condition is described by the Euler equations. This suggests 
that, unless X(u) —y 0, the Euler equations (8) might also fail to govern 
the inviscid limit of the Navier-Stokes equations with slip boundary 
conditions (1). 

3. From the Boltzmann Equation with Accomodation 
Boundary Condition to the Navier-Stokes Equations 
with Slip Boundary Conditions 

In this section, we revisit the incompressible Navier-Stokes limit for 
the Boltzmann equation in the case of the initial-boundary value prob- 
lem. Our main interest is to understand how the slip boundary condi- 
tion arises from Maxwell's accommodation boundary condition at the 
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kinetic level in the fluid dynamic limit, and especially how the slip co- 
efficient is related to the accommodation parameter. Strictly speaking, 
this is not needed in the proof of the main result in the present paper. 
Therefore, the discussion in this section will be only formal, along the 
line of [4]. 

Consider the Boltzmann equation with the incompressible Navier- 
Stokes scaling 

(15) ed t F e + v-V x F e = -B(F ei F e ). 

Here the unknown is the distribution function F = F(t, x, v) that is 
the density at time t of molecules with velocity v G at the position 
x G Q with respect to the phase space Lebesgue measure dxdv. 



3.1. Formal structure of the Boltzmann equation. The Boltz- 
mann collision integral acts only on the v variable in F e , keeping t, x 
as parameters. Its expression for <p G C C (R N ) is 
(16) 

B((j>, <f>)(v) — // (<f>(v')(f>(vl) — (f>(v)(f>(v*))b(v — v*,u)dv*dw 

where v', v'^ G are the velocities of 2 identical particles about to 
undergo an elastic collision, assuming that their post-collision velocities 
are v,v* G H N . The set of all possible pre-collision velocities v\v[ are 
parametrized by the unit vector co as follows: 

v = v'(v, v*, to) :— v — (v — v*) ■ coco , 
v'„ = v'„(v : w*, to) := v* + (v — v*) ■ coco . 

The collision kernel b(z, to) > is a locally integrable function that 
satisfies the symmetries 

(18) b(v — i;*, to) = b(v* — v, to) = b(v' — v'„ : to) 

a.e. in (v,v*,u), assuming that v' and v[ are given in terms of v,v*,co 
by the relations (17). Depending on the growth of the collision kernel b 
as |i> — v*\ — > +oo, the collision integral can be extended by continuity 
to larger classes of functions than C C (R N ). Finally, we denote 

B(F, F)(t, x, v) := B(F(t, x, •), F(t, x, -))(v) . 



(17) 
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The collision integral satisfies the identities 



(19) 



B((f),(f))(v)dv = 0, 
B((f),(f))(v)vdv = 0, 



B(<p,<i))(v)\v\ 2 dv = 0, 



for each G C^R^) or in the larger class allowed by the growth at 
infinity of the collision kernel b. As a result, whenever F is a classical 
solution of the scaled Boltzmann equation (15) with appropriate decay 
as \v | — > oo, 



(20) 



ed t / + div-r / vF e dv = , 

Jn N Jr n 

ed t [ vFJv + div x [ v m F e dv = , 
ed t / \\v\ 2 F t dv + div x / v\\v\ 2 F e dv = , 

V JRN JrN 



and these relations are the local conservation laws of mass, momentum 
and energy respectively. 

Whenever b{v — v*,ou) has polynomial growth as \v — v*\ — > oo, for 
each positive, rapidly decaying G C(R JV ) such that ln0 has polyno- 
mial growth as — > oo, Boltzmann's H Theorem states that 



(21) 
and 

(22) 



B((f),(f))(v)\n(f)(v)dv < 0, 



£(0, 0) (u) In 0(y)cfo = <£> B{(j>, 0) = 



R J 



is a Maxwellian distribution, 
meaning that there exists p, 9 > and u G R w such that 



(23) 



<j>{v) = M p , u ,e(v) := 



{2n0) N / 2 



e 26 



for all v G R^. As a result, whenever F is a classical solution of the 
scaled Boltzmann equation (15) with appropriate decay as \v \ — > oo, it 
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satisfies the differential entropy inequality 

ed t / F e In F e dv + div x / vF e \nF e dv 

(24) J * N 

= — / B(F e ,F e )\nF e dv < 0. 
e Jr n 

Throughout this paper, we denote 

(25) M:=M lfi ,i. 

Since Maxwellians are equilibrium distributions for the collision in- 
tegral, it is natural to investigate the linearization thereof about a 
Maxwellian, say M for simplicity — the case of an arbitrary Maxwellian 
being similar. We therefore introduce the linearized collision operator 
in the form 

i.e. 

£m4>{v) '■= // {4>{v) + 4>{v*)— <p{v') — (p(v'J)b(v— v*,u)Mdv*du . 

Under certain assumptions on the collision kernel b, known as Grad's 
angular cutoff assumption, H. Grad proved in [17] that Cm is an un- 
bounded, self-adjoint Fredholm operator on L 2 (R N ; Mdv) with domain 

D(£ M ) := {0 G L 2 (R N ; Mdv) | <f>(b* v M) G L 2 {R N ; Mdv)} , 

where 

b{z) := / b(z, u)du 



and * v designates the convolution product in the v variable. Moreover, 
the nullspace of Cm is 

Ker Cm = span{l, vi, . . . , v n , \v\ 2 } . 

In particular, the tensor field A{v) = v® 2 — jy\v\ 2 satisfies A±Ker Cm, 
so that, by the Fredholm alternative, there exists a unique tensor field 

A G D(£ M ) n (Ker C) L such that C M A = A 

componentwise. 

Henceforth in this paper, we assume that the collision kernel b comes 
from a hard cutoff potential in the sense of Grad, and more precisely 
that it satisfies, for some constant C b > and all (z,u) G x S^ -1 , 

(26) < b{z, uj) < C b (l + \z\), and b(z) > . 
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3.2. Boundary value problem and fluid dynamic limit. The in- 
compressible Navier-Stokes limit of the Boltzmann equation bears on 
solutions of the Boltzmann equation that are of the form 

(27) F e = M(l + eg e ), 

where it is understood that the relative number density fluctuation g t 
is 0(1) in some sense to be made precise as e — > 0: see [3, 4, 5] for more 
details, together with physical justifications for this scaling assumption. 

Here, the scaled Boltzmann equation (15) is set on the spatial domain 
Q, with Maxwell's accommodation at the boundary, that is assumed to 
be maintained at the constant temperature 1. This boundary condition 
reads 

F £ (t, x, v) = (l- a)1Z x F e (t, x, v) + aA x (J^J (t, x)M{v) , 

x G dfl , v ■ n x < , 

where 

(29) 7Z x F(t, x, v) := F(t, x, v — 2v ■ n x n x ) 

is the specular reflection operator on the boundary, while 

(30) A x cf):=V2^[ cf)(v)(vn x ) + M(v)dv. 

Jr n 

In (28), the parameter a satisfies < a < 1, and is called the ac- 
commodation coefficient. The case a = corresponds with specular 
reflection of the gas molecules on dfl without thermal exchange, while 
the case a = 1 corresponds with diffuse reflection, or total accommo- 
dation, in which case gas molecules are instantaneously thermalized at 
the boundary. 

Henceforth, we denote, for each e L 1 (R N ; Mdv), 



(0) := / ip(v)M(v)dv. 
Jr n 



Theorem 3.1. Let (F e ) t>0 be a family of solutions of the scaled Boltz- 
mann equation (15) set on the spatial domain Q, satisfying the accom- 
modation boundary condition (28) on dfl. Assume that the relative 
fluctuations 

F e — M 

a.e. and in weak-Lj oc (Il + xQx~R N ; dtdxMdv) (possibly up to extraction 
of a subsequence), and that 

(\v\H\ v \ >R \g t \) + (\A\l lvl>R \Q(g e , g e )\) + (|i|Ml H>i? |^|> -> 
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in Lj oc (R + xO) as R — >■ +00 uniformly in e > 0. T/ien 

f ) = p(£, x) + x) • v + 9(t, x)^(\v\ 2 — N) , 
where u is a solution of the incompressible Navier-Stokes equations 
div x u = , 
^ d t u + div^w® 2 ) + V x p = vA x u , 
and where 

V = (N-l)(N+2) (^ '■ A) . 

Assume further that 

(v-n x g £ \ dn ) ->• (u-n^lan) 
(u T (u • n x ) + g e \ m ) ->■ (u T (u • n^^) 

in weaA;-L z 1 oc (R + x cftl) . T/ien, i/ie velocity field u satisfies the boundary 
condition 

u ■ n x = , x G <9f2 , 



z/(£(n) • ra x ) T + Am = , x G <9f2 , 
where the slip coefficient is given by the formula 



2tt 

Proof. Observe that, under the substitution i> 1— >• w = t> — 2t> • nn, one 
has, for each unit vector n, 



4>(v — 2v ■ nn)(v ■ n)-Mdv = / <fi(w)(w ■ n) + Mdw , 

Jr^ 

so that, for each x G <9f2, one has 

/ F e (t, x,v)v ■ n x dv = —a / F £ (t, x, v )(v ■ n x ) + dv 
Jn, N Jr. n 

+ aA x (J^J (t,x) J ^M(v)(vn x )-dv = 0, 

since 

/ M(v)(v ■ n x ) + dv = I M(v)(v -n x )-dv = — ]= . 

JR'V JR^ V27T 

Hence, for each x G <9f2 and each e > 0, one has 

(vg t )(t,x) ■ n x = - / F e (t, x, v)v ■ n x dv = , 

so that, after passing to the limit as e — > 0, 

u(t, x) ■ n x — (vg)(t, x) • n x — , £ > , x G <9f2 . 
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Next write the local conservation of momentum — the second local 
conservation law in (20) in the form 

d t (vg e ) + div x ^(Ag e ) + V x -^\v\ 2 g e ) = , 

where 

A = A(v) := v® 2 - ±\v\ 2 . 

Let now w = w(x) G H N designate a compactly supported C 1 vector 
field on H N satisfying 

divw = , and w(x) ■ n x — , x G dfl . 

Taking the inner product of both sides of the local conservation of 
momentum with w and integrating over fl leads to 

d t I w ■ (vg e )dx + / w <S> n x : -(Ag e )dS(x) 
Jn Jdn e 

- ! Vw : -(Ag e )dx = 0, 
Jn e 

J w -V x ^(jf\v\ 2 g e )dx = J div x (^-(±\v\ 2 g e ) ■ w^j dx 

= / -(jf\v\ 2 ge)w-n x dS{x) = 0. 
Jdn e 

Next we pass to the limit in each term appearing in (31): following 
the analysis in [3, 4], one finds that 

/ w ■ {vg € )dx — > / w • (vg)dx = w ■ udx 
Jn Jn Jn 

while 

/ Vw : -(Ag e )dx 
Jn e Jn 



(31) 



since 



Vw : {u m - vll(u))dx. 

(Indeed, since w is divergence-free, Vw : ^|m| 2 = ^|u| 2 div x u = 0.) 
It remains to analyze the boundary term 

/ w ® n x : -(Ag e )dS(x) . 
Jan e 

Since w ■ n x = on dQ, 

w®n x : (Ag e ) = w ® n x : (v m g e ) = (v T v ■ n x g e ) ■ w . 
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At this point, we decompose the boundary term into the contribution 
of gas molecules about to collide and those having just collided with 
the boundary 

(v T v ■ n x g e ) = (v T v ■ n x l v . nx>0 g e ) + (v 

and use the accommodation condition to write 

{v T v ■ n x l v . nx<0 g e ) = (v T v ■ n x l v . nx<0 ((l - a)TZ x g e + aA x (g e ))) . 

Observing that 

(v T v ■ n x l v . nx<0 TZ x g e ) = (TZ x (v T v ■ n x l v . nx>0 )g e ) 

= -{v T v ■ na.l w . nil ,>o((l - a)g e + aA x (g e ))) , 

we conclude that 

(v T v ■ n x l v . nx<0 g e ) = -(v T v ■ n x l v . nx>0 ((l - a)g e + aA(g e ))) 
so that 

{v T v ■ n x g e ) = a(v T v ■ n x l v . nx>0 (g e - A x (g e ))} = a(v T v ■ n x l v . nx>0 g e ) 

- where the second equality follows from the fact that the function 
v^v T v- n x l v . nx>0 (g e - A x (g e ) is odd in v T . 

Therefore the boundary term appearing in (31) becomes 

/ w <S> n x : -(Ag e )dS(x) = — f w ■ (v T v ■ n x l v . nx>0 g t )dS(x) . 
Jan e e Jan 

Assume that a = a(e) varies with e so that a(e)/e — > a as e — >■ 0. 
Since 

g e -+ g = p + u ■ v + 6\(\v\ 2 - N) 
and we already know that 

u • n x = on dfl , 

one has 

w ■ {v T v ■ n x l v . nx>0 g e ) -> (vf 2 {v ■ n x )+) :u t ®w 

= jtt(\vt\ 2 {v -n x )+)u-w . 
Thus, passing to the limit in (31) leads to 



(32) 



d t / w ■ udx + jfzi(\v T \ 2 {v ■ n x ) + ) I u-wdS(x) 
Jn '' Jan 

- I Vw : {u® 2 - vY,{u))dx = . 
Jn 



(Notice that the term (|ti T | 2 (ti -n x ) + ) is independent of x and therefore 
comes out of the boundary integral.) 
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Whenever u(t, •) G C 2 (Q), applying Green's formula transforms the 
last integral above into 

/ Vw : (u® 2 - vY,{u))dx = - I W (div x (u m ) - vdiv x (Z(u)))dx 
Jn Jn 



+ I w • uu ■ n x dS{x) 
'an 



— v \ w ■ • n)dS(x) . 

JdQ 

Since u-n x = on dQ, the second integral on the right-hand side above 
vanishes, and since div^-u = 0, one has div x (T,(u)) = A x u, so that 

/ Vu> : {u® 2 — z/£(w))<ir = — w ■ (div x (u® 2 ) — uA x u)dx 
Jn Jn 

-v w (E(u) • n)dS(x) . 
Jon 

Thus, if u G C 2 ([0,T] x H), the equality (32) becomes 

d t w ■ udx + w ■ (div x (u® 2 ) — uA x u)dx 
Jn Jn 

+ jTi(\vt\ 2 {v -n x ) + ) I u-wdS(x)+v I w(E(u)-n)dS(x) = . 



ion Jan 

This identity holds, say, for each w G C£°(fi;R iV ). In particular, it 
holds for each w G C£°(Q; H N ), which implies that 

d t u + dw x (u® 2 ) - vA x u = -V x p 

in the sense of distributions (for some p G X>'(R+ x ^))- Since the 
velocity field u G C 2 ([0,T] x Tt) we conclude that p e_C\[0,T] x H). 
Substituting this in the identity above with w G C™(Q; ~R N ) gives 

— / w ■ V x pdx + / (z/£(w) • n x + jfzi(\v T \ 2 {v ■ n x ) + )u) ■ wdS(x) = 
Jn Jan 

and since, by Green's formula, 

— / w ■ V x pdx = — / div x (pw) dx = / pw ■ n x dS(x) = , 
Jn Jn J an 

we conclude that 

i/(E(u) • ra x ) T + £o-(\ Vt \ 2 (v ■ n x ) + )u = onffl. 
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In other words, (32) is the weak formulation of 

d t u + div x (u® 2 ) - uA x u = -V x p, x e Vl , t > , 
v(E(u) ■ n x ) T + Am = , x G dVL , t>0, 

u • n x — , x G <9f2 , t > . 

with 

A = # a r(HVn,) + > ° n 



2tt 

□ 

The argument above is a proof of the (formal) Navier-Stokes limit 
Theorem 3.1 by a moment method analogous to the one used in [4]. 
As far as we know, the first derivation of this slip boundary condition, 
in the steady, linearized regime (i.e. leading to the Stokes equations 
in the fluid limit), is due to K. Aoki, T. Inamuro and Y. Onishi [1] 
(see especially formula (33) in that reference). That derivation uses 
a Hilbert expansion method (formal series expansion of the solution 
of the Boltzmann equation in powers of the Knudsen number e). The 
interested reader is referred to the recent book by Y. Sone [35] (in 
particular to §3.7 there) for a systematic study of boundary conditions 
in the context of the fluid dynamic limit of the Boltzmann equation. 

For a complete proof of the derivation of the same slip boundary con- 
dition as in Theorem 3.1 in the linearized regime — i.e. in a situation 
where the limiting equation is the Stokes, instead of the Navier-Stokes 
equations — the reader is referred to the work of N. Masmoudi and L. 
Saint-Raymond [26]. 



4. From the Boltzmann Equation with Accomodation 
Boundary Condition to the Incompressible Euler 

Equations 

In this section, we consider the Boltzmann equation in the scal- 
ing leading to the incompressible Euler equations in the fluid dynamic 
limit. We recall from [3, 4] that this scaling is 

(33) ed t F e + v ■ V x F e = -J^B(F e , F e ) , (x,v)eQxR N , 

with q > 0, while the distribution function F e is sought in the same form 
(27) as in the Navier-Stokes limit. This scaled Boltzmann equation is 
supplemented with Maxwell's accommodation condition on dfl, with 
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accommodation coefficient a = a(e) driven by the small parameter e: 

F e (t, x, v) = (l- a(e))K x F e (t, x, v) + a(e)A {^j (t, x)M{v) , 

x G dQ , v • n x < , 

and with the initial condition 

(35) F £ (0,x,v) = F™(x,v), (x, v) E Q x ~R N . 

The formal result presented in Theorem 3.1 suggests that, in the 
limit as e — > 0, the velocity field 

lim - / vF f dv 



e Jr n 

should behave like the solution of the incompressible Navier-Stokes 
equations with kinematic viscosity of order e q and with slip bound- 
ary condition with slip coefficient of the order of a(e)/e. Thus, if 
a(e) = o(e), Theorem 2.2 suggests that this velocity field should satisfy 
the incompressible Euler equations (8). In fact, the formal result in 
Theorem 3.1 is only a guide for our intuition, and we shall give a direct 
proof of the Euler limit starting from the Boltzmann equation with 
accommodation boundary condition without using the Navier-Stokes 
limit. 

4.1. Renormalized solutions and a priori estimates. Global so- 
lutions of the Cauchy problem for the Boltzmann equation for initial 
data of arbitrary size have been constructed by R. DiPerna and P.-L. 
Lions [12]. Their theory of renormalized solutions was extended to 
the initial boundary value problem by S. Mischler [27]. His result is 
summarized below — see also section 2 in [26] and section 2.3.2 of [31]. 

Theorem 4.1 (Mischler). Let F™ = F™(x,v) > a.e. on Q x ~R N be 

a measurable function satisfying 

(1 + M 2 + \\nFi n (x,v)\)Fi n (x,v)dxdv < +oo . 

QxR N 

There exists F t e C(R + ;L 1 (fi x H N )) satisfying the initial condi- 
tion (35), and the Boltzmann equation (33) together with the bound- 
ary condition (34) in the renormalized sense, meaning that, for each 
T G C^R-j-) such that Z h-> \f\ + ZT'(Z) is bounded on R +; the func- 
tion 

r ' (S) m > Fe) e LL{R+ X ^ X RN) 
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and 

III r ( S J ( e(9 * + u ' V x )<pMdvdxdt 
= L II 1 ^ ' nxMdvdS ^ dt 

- t /L. r ®*i« M * ,fa 

/or eac/i G C C 1 (R + x x R^). 

Moreover 

t/ie irace of F t on dfl satisfies the accommodation boundary condition 

F e \ dn (t,x,v) = (1 - a)1Z x (F e \ m )(t, x, v) + A x (^-jf^J (t,x)M(v) 

for a.e. (t, x, v) G R+ x dQ x R N such that v ■ n x > 0; 
b) the distribution function F t satisfies the local conservation law of 
mass 



ed t / F e dv + div^ / vF e dv = 
Jn N Jn N 

with boundary condition 



/ F e (t,x,v)v 
Jn N 



n x dv = , x G dQ , t > ; 



c,) t/ie distribution function F e satisfies the relative entropy inequality 
H(F e \M)(t) - H(F™\M) 



< 



— — I I V t (s,x)dxds / I VQ e (s,x)dxds 

e 9 Jo Jn e Jo Jan 



for each t > 0, where the following notations have been used: for each 
f, g measurable onQx H N such that f > and g > a.e. , the relative 
entropy is 

H(f\g) := ff h(--l) 9dvdx 

with 

h(z) := (1 + z) ln(l + z) - z , 
while the entropy production rate per unit volume is 

V e := HI r ( - 1 J F e F e *b{v - v*,u)dvdvju 



22 C. BARDOS, F. GOLSE, AND L. PAILLARD 

with 

r(z) := z\n(l + z) > , 
and the Darrozes-Guiraud information is 



In particular H(f\g) > since h > on [— 1, +oo) and V e > a.e. on 
R + x Q since r > on (— l,+oo), while T>Q e > a.e. on R + x dfl 
since h is convex and A x is the average with respect to a probability 
measure; 

d) for each T > and each compact K C dVt, there exists Ck,t > 
such that, for each e > 0, one has 

I II F e (t,x,v)(v-n x ) 2 M(v)dvdS(x)dt<C KtT , e>0. 

JO JjKxR N 

Statement d) appears in [26], without proof. We give a brief justifi- 
cation for this estimate below. 

Notice that in general, renormalized solutions of the initial-boundary 
value problem (33)-(34)-(35) are not known to satisfy the local conser- 
vation law of momentum — see equation (2.35) in [31] for a variant 
involving a defect measure, following an earlier remark due to P.-L. 
Lions and N. Masmoudi. 

At variance, any classical solution F e G C(R + x Q x H N ) is contin- 
uously differentiable in (t, x) and such that 

v ^ sup (\F e (t,x,v)\ + \d t F e {t,x,v)\ + \V x F e {t,x,v)\) 

, , 0<t<T 

(36) \*\<R 

is rapidly decaying as |i>| — >■ +oo , 

one has 

ed t / vFJv + div x / v m F € dv = , x e £1 , t > . 
Jr n Jr n 

Moreover, for each w G C^(R + x fi), Green's formula implies that 

/ / / (ev • d t w(s, x) + v® 2 : V x w(s, x))F e (s, x, v)dvdxds 
Jo JJnxn N 



(37) 



J JJixR" 

/ v® 2 : w(s, x) ® n x F e (s, x, v)dvdS(x)ds 

JdttxR N 

+e / / w(t, x) ■ vF e (t, x,v)dvdx 
J JnxR N 



I In 



w(0, x) ■ vF* n (x, v)dvdx . 

nxn N 
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Let us use the accommodation condition (34) to reduce the boundary 
integral: 

/ v ® 2 : w(t, x) Cg) n x F e (s, x, v)dv 
w(s, x)-v)(v-n x ) + F e (s, x, v)dv 

R ' 

[w(s, x)-v)(v-n x )-F e (s, x, v)dv 
and, whenever w is tangential on dfl, one has 

(w(s, x) -n x )~F e (s, x, v)dv 
= J ^(w(s,x)-v)(v-n x )- ^(1 - a)K x F e + aA x (s,x,v)dv 
= J ^(w(s,x)-v)(v-n x ) + ^(1 -a)F t + aA x (J^J M^j (s,x,v)dv. 

Therefore 

/ v® 2 : w(s, x) <g> n x F e (s, x, v)dv 
Jn N 

— a (w(s, x) ■ v )(v ■ n x ) + F e (s, x, v )dv 
-«Ax {~^J J n (w(s,x) ■ v)(v ■ n x )+Mdv 

and the last integral vanishes since the integrand is odd in the tangen- 
tial component of v. 

Finally, whenever w G C*(R+ x Q) is tangential on dQ, one has 

/ v® 2 : w(s, x) ® n x F e (s, x, v)dv 
Jr n 

— a (w(s, x) ■ v )(v ■ n x ) + F e (s, x, v )dv . 
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Therefore, each classical solution F t G C(R + x Q x R w ) of the initial- 
boundary value problem (33)- (34)- (35) that is continuously differen- 
tiable in (t, x) and satisfies (36) also verifies 



Henceforth, we shall consider exclusivey renormalized solutions of the 
initial-boundary value problem (33)- (34)- (35) satisfying the identity 
(38) for each vector field w G C£(R+ x fi) tangential on <9f2 and such 
that divj, w = on Q. 

Now for estimate d) in Theorem 4.1. 

Proof of estimate d). A renormalized solution of the initial-boundary 
value problem (33)-(34)-(35) can be constructed as the limit for e > 
fixed and m — > +oo, of solutions F e m of the approximating equation 



with the same initial and boundary conditions (34)- (35) satisfied by 
F e m , where the approximate collision integral is given by the same 
expression as Boltzmann's collision integral with collision kernel b re- 
placed with its truncated variant b m defined as 1 



Let £ be a compactly supported C 1 vector field satisfying 

£(x) = a(x)n x for each x G dVL , with a > on dfl and a = 1 on K . 

Since the approximate collision integral in (39) is normalized with an 
average of F e ^ m with respect to v, all solutions of that equation satisfy 
equality (37) for any w G C^(R + x R^), i.e. the local conservation of 



1 For all x, y € R, the notation x Ay designates min(x, y). 




: V x w(s,x))F e (s,x,v)dvdxds 



(38) 





b m (v — v*, of) := m A b(v — v*, u) . 
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momentum. In other words, 

III (v m : V^(:r))F e , m ( S , x, v)dvdxds 
Jo JJqxr n 

= // (v-Z(x))(v-n x )F e>m (s,x,v)dvdS(x)ds 

J0 JJdflxR N 

+e / / • vF e>m {t, x, v)dvdx 

J JnxR N 

-e £(x) ■ vF™ m (x, v)dvdx . 

J JnxR, N 



Therefore, since F em > a.e. and (v ■ £(x))(v ■ n x ) = a(x)(v ■ n x ) 2 > 
for each x G dQ and v G R^, one has 



0< / // (v ■ n x ) 2 F €im (s, x,v)dvdS(x)ds 
Jo JJkxr n 

/ // (v ■ £(x))(v •n x )F ejm (s,x,v)dvdS(x)ds 
Jo JJdnx-R N 

= 111 (*>® 2 : Vxt(x))F etm (s, x, v)dvdxds 
Jo J Jnxn N 



< 



(40) 



+6 // ^x)-vF^ m (x,v)dvdx 

J JnxR N 

1 / / €( x ) ' v F,m(t, x, v)dvdx . 

J JnxR N 



At this point, we recall that the function h : z \— > (1 -\- z)\n{l -\- z) — z 
introduced in Theorem 4.1 has Legendre dual 



h*({)= sup((z-h(z)) = e<-(-l. 

Z>-1 

By Young's inequality — or equivalently, by definition of h* (see for 
instance [5]) — one has 

i(l + \v\ 2 )e\g e , n \ < h(e\g e>n \) + h*(\{l + \v\ 2 )) 
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so that, for each nonnegative x £ C C (FL N ), 

/ / x(x)(l + \v\ 2 )F e)m (t,x,v)dvdx 
J Jnxn N 



[ [ x(x){l + M 2 )(l + eg e>m )(t,x,v)M(v)dvdx 

J JflxR N 



(41) < // x(x)(l + h*(\(l + \v\ 2 ))M(v)dvdx 

JjQxR N 

+ \\ X \\L<~H(F €;m \M)(t) 



< // x(x)(l + h*(\(l + \v\ 2 ))M(v)dvdx 

+ \\x\\ L °°H(F: n \M) . 

since the relative entropy estimate c) in Theorem 4.1 is also satisfied 
by the approximate solution F e m . 

Since £ is compactly supported, putting together (40), (41), and 
letting m — > +oo leads to estimate d). □ 

4.2. The Euler limit. Let u in € H{$1), and pick initial data F™ for 
the Boltzmann equation satisfying 

(42) ^H(F?\M lt€ui n tl ) 

Theorem 4.2. For eac/i e > 0, Zei F e &e a renormalized solution of 

(33) -(34)-(35) satisfying the local momentum conservation law (38) for 
each w E C^(R + x Q) satisfying div x w = and w -n x \ m = 0. Assume 
that the accommodation parameter a in the accommodation condition 

(34) at the boundary depends on the scaling parameter e in such a way 
that 

a(e) = o(e) as e — > . 
Then, for each compact K C Vt, the family 

1 f 

- / vF e dv 
e Jr n 

is relatively compact in L°°(R + ; ^(K)), and each of its limit points as 
e — > is a dissipative solution of the Euler equations (8). 

Assume that u m G "H(f2) is smooth enough so that the initial- 
boundary value problem for the Euler equations (8) has a classical 
solution on some finite time interval [0, T] — for instance u m G H S (Q) 
with s > y + 1, or u m G C 1 ' 9 with < 9 < 1. In that case, the con- 
vergence result above can be strengthened with the notion of entropic 
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convergence, invented by Dave Levermore specifically to handle such 
problems. 

Definition 4.3 (C. Bardos, F. Golse, CD. Levermore [5]). A family 
g e = g e (x,v) of Lj oc (Q x H N ; Mdxdv) is said to converge entropically 
at order e to g = g(x,v) as e — > if the following conditions hold 

(i) 1 + eg e > a.e. on Q x R w for each e, 

(ii) g e —> g weakly in L} oc (Q x H N ; Mdxdv) as e — > 0, 
(Hi) and 

\H(M(l + eg e )\M) ->■ \ [[ g(x,v) 2 M(v)dxdv 
e JJqxr n 

as e — > 0. 

We recall that, if g e — )■ g at order e, then, for each compact K C fi, 
one has 

(1 + \v\ 2 )\g e (x,v) - g(x,v)\M(v)dvdx -^0 as e ->■ . 

'kJb, n 

In other words, entropic convergence implies strong L 1 convergence 
with the weight (1 + \v \ 2 )M(v) (see Proposition 4.11 in [5].) 

Whenever the incompressible Euler equations (8) have a classical 
solution u on [0, T] x Q, using the weak-strong uniqueness property of 
dissipative solutions and the conservation of energy satisfied by classical 
solutions of (8), we arrive at the following stronger convergence result, 
which is a straightforward consequence of Theorem 4.2. The interested 
reader is referred to the proof of Theorem 6.2 in [5] where the squeezing 
argument leading from weak compactness to entropic convergence is 
explained in detail. 

Corollary 4.4. Consider a family F™ = F™(x,v) > a.e. of mea- 
surable functions on Q x H N such that 

Fi n {x,v)-M{v) in( . 

. . — — -> u m (x) ■ V 
eM(v) 

entropically of order e as e — > 0, where u m G Hi^t) is smooth enough 
so that the initial-boundary value problem (8) has a classical solution 
u defined on the time interval [0,T] with T > 0. 

For each e > 0, let F e be a renormalized solution of (33)-(34)- 
(35) satisfying the local momentum conservation law (38) for each 
w E C^(R+ x Q) that satisfies div x w = and w ■ n x \ m = 0. Assume 
that the accommodation parameter a in the accommodation condition 
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(34) at the boundary depends on the scaling parameter e in such a way 
that 

a(e) = o(e) as e — > . 

Then 

F e (t,x,v)-M(v) 

— — > u(t,x) ■ V 

eM(v) v 1 

entropically of order e as e — > 0, for a.e. t e [0,T]. 

The proof of Theorem 4.2 above occupies the remaining part of the 
present section. 

4.3. The relative entropy inequality. Statement c) in Theorem 4.1 
bears on the evolution of the relative entropy of the distribution F e 
with respect to the uniform Maxwellian M — A4 1,0,1 ■ In the next 
proposition, we consider the evolution of the relative entropy of the 
distribution F e with respect to a local Maxwellian of the form A^i, eW) i, 
where w is a solenoidal velocity field on il that is tangential to dil. 

Proposition 4.5. Let w G C C 1 (R + x Q) be such that 
div x w = and w ■ n\ ^ = . 



an 



Then, for each e > 0, renormalized solution F t of the initial-boundary 
value problem (33)-(34)-(35) satisfying the momentum conservation 
identity (38), also satisy the relative entropy inequality 
(43) 

j 2 H{F e \M w ){t) - ^H(F^ n \Mi : ew(o,-),i) 

< — t— / / V e (s,x)dxds / / VQ e (s,x)dxds 

e Jo Jn e Jo Jan 

/ II (v — ew(s,x))® 2 :V ' x w(s,x)F e (s,x,v)dxdvds 

e Jo JJnxR N 

/ // (v — ew(s, x)) ■ E(w)(s, x)F e (s, x, v)dvdxds 

e Jo J JnxR N 



a 



+ — / / (w(s, x) ■ v)(v ■ n x ) + F e (s, x, v)dvdS(x)ds 



t 2 Jo Janxn N 



for each t > 0. 
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Proof. We begin with the straightforward identity 
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H{F e \M w ) = H(F e \M) + [[ FJn ( T ^—] dxdv 

JJnxR, N VAli, eti ,,i/ 

= H(F t \M)+ [[ F e (\\v - ew\ 2 - \\v\ 2 )dxdv 
JJnx~R N 

= H(F e \M)+ [[ F e (\e 2 \w\ 2 - ev ■ w)dxdv . 
J JnxR N 



Thus 

(44) 

H{F e \M w )(t) - H{F e \M w ){Q) = H(F e \M)(t) - H(F e \M)(0) 

+\e 2 / / F e (t, x,v)\w(t, x)\ 2 dxdv 

-\e 2 ff F* n (x,v)\w(0,x)\ 2 dxdv 
J Jnxn N 

— e / / F e {t, x,v)v ■ w(t, x)dxdv 
J Jnx-R N 



+e 



/ / F™(x, v)v ■ w(0, x)dxdv . 

J JnxR N 



According to the continuity equation in statement a) in Theorem 4.1 



+ \e 2 



// F e (t, x, v )\w(t, x)\ 2 dxdv 
JJnxn, N 



-\e 2 / / F* n (x,v)\w(0,x)\ 2 dxdv 

, JJnxR N 

(45) ft fj . 

= // F e (e 2 d t + ev ■ V x )\\w\ 2 dxdvds 
Jo J Jnxn N 

— / / F e w ■ (e 2 d t w + ev ■ V x w)dxdvds . 
Jo J Jnxn N 



In (44), we replace the term H(F e \M)(t) - H(F e \M)(0) with the right 
hand side of the inequality of statement c) of Theorem 4.1, and use 
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(45) together with (38) to arrive at 
(46) 

H(F e \M 1 , ew , 1 )(t) - H(F e \M 1 , ew , 1 )(0) 

< — ^— J I V e (s,x)dxds I J VQ e (s,x)dxds 

e q Jo Jn e Jo Jan 



ft 

+ 



/ // w ■ (e 2 d t w + ev ■ V x w)(s, x)F e (s, x, v )dxdv ds 
Jo J JnxR N 

— // (ev ■ d t w(s, x) + v® 2 : V x w(s, x))F e (s, x, v)dvdxds 

Jo J Jnx~R N 

+ a / (w(s, x) ■ v)(v ■ n x ) + F e (s, x,v)dvdS(x)ds . 

Jo JdQxR N 

Next we express d t w in terms of E(w) = d t w + w ■ V x w and V x w: thus 

(ev ■ d t w(s, x) + v ® 2 : V x w(s, x)) — w ■ (e 2 d t w + ev ■ V x w) 
= (v - ew)® 2 : V x w + e(v - ew) ■ E(w) . 

In the right hand side of (46), we substitute 

/ // w ■ (e 2 d t w + ev ■ V x w)(s, x)F e (s, x, v)dxdvds 
Jo JJq.xr n 

— // (ev ■ d t w(s, x) + v 02 : V x w(s ) x))F e (s, x, v )dv dxds 
Jo J JnxR N 

= — // (v — ew(s, x))® 2 : V x w(s, x)F t (s, x, v )dxdvds 

Jo J JnxR N 

el II (v — ew(s, x)) ■ E(w)(s, x)F e (s, x,v)dvdxds 
Jo J JnxR N 

and arrive at the relative entropy inequality (43). □ 

4.4. Control of the boundary term. The relative entropy inequal- 
ity (43) is the same as in the one considered in [10], [25] and [29], 
except for the boundary term — i.e. the last term on the right hand 
side, which is in general not nonpositive. Since the effect of the bound- 
ary is our main interest in this paper, and the only difference with the 
case of the Cauchy problem treated in [29], the core of our argument 
is to obtain a control of that term. 

Lemma 4.6. With the notations of Theorem for each e > and 
each w G C*(R + x Q) that is tangential on the boundary dfl, each 
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renormalized solution F t of the initial-boundary value problem (33)- 
(34)-(35) satisfies the inequality 



a 

72 



(w(s, x) ■ v)(v • n x ) + F e (s, x, v)dv 



(47) =— (w(s,x) ■ v)(v ■ n x ) + g e (s,x,v)M(v)dv 

e Jn N 

ol ol 
< -^VQ e (s,x) + -C(w)l w{s/x)f L A x (F e )(s,x) 

a.e. in (s,x) G R + x dQ, where 

C(w) := \ I ( e 2 IMU~M - 2||iu|| L oo|u| - l)( Vl )+M(v)dv . 
Jr n 

We use Young's inequality for a translate of the function h defined 
in Theorem 4.1, much in the same way as in the proof of Theorem 6.2 
in [5] (see especially pp. 738-739 there). 

Proof. Let z > — 1; for each z > — 1, set 

l(z - z Q ) := h(z) - h(z ) - h'(z )(z - z ) ■ 

We recall that the Legendre dual of the function h defined in Theorem 
4.1 is 

h*(p) := sup (pz — h(z)) = e p — p — 1 , pGR. 

z>-\ 

A straightforward computation shows that 
l*{p) = sup (p(z - z ) - l(z - z )) 

Z>-1 

= sup (p(z - z Q ) - h(z) + h(zo) + h'(z Q )(z - z Q )) 

z>-l 

=h(z ) - (h'(z ) + p)z + sup ((h'(zo) + p)z - h(z)) 

Z>-1 

=h(z ) - (h'(z )+p)z + h*(h'(z )+p) 

=h(z ) - (h'(z ) + p)z + e h '^e p - h'(z ) -p-1 
= (1 + z ) ln(l + z ) - z - (ln(l + 2 ) + p)^ 

+ (1 + 2 )e p - ln(l + 20) - P ~ 1 
= (l + ^ )(e p -p-l) = (l + ^o)^(p). 
Writing F e = M(l + eg e ) and observing that 

/ (w(s,x) ■ v)(v ■ n x ) + M(v)dv = , (s, x) G R+ x <9fi 



32 C. BARDOS, F. GOLSE, AND L. PAILLARD 

since w is tangential on dQ, one has 



(48) 



/ (w(s, x) ■ v)(v ■ n x ) + F € (s, x,v)dv 

(w(s, x) ■ v)(v ■ n x ) + g e (s, x, v)M(v)dv 



Then 

/ (w(s,x) ■ v)(v ■ n x ) + g e (s,x,v)M(v)dv 
Jr n 

= / (w(s,x) ■ v)(v ■ n x ) + (g e (s,x,v) - A x (g e ))M(v)dv 
Jr n 

by (48) since A x (g e ) is independent of v. 

By definition Young's inequality (or equivalently, by definition of the 
Legendre dual of I), 

2e 2 (g e - A x (g e ))(w ■ v) < l(e(g e - A x (g e ))) + l*(e(2w ■ v)) , 

so that 

2 1 (wv)(v n x ) + g t Mdv < \ [ l{e{g e - A x (e)))(v ■ n x ) + Mdv 

JR N e JR N 

+— / l*(e{2w v)){v ■ n x ) + Mdv . 



e 2 ./ R iv 



First, since A x is the average under a probability measure, one has 
l(e(g e -A x (e)))(vn x )+Mdv 



IR N 

= ^A x (l(e(g e - A x (e)))) 

= -±=A x (h(eg e ) - h(eA x (g £ )) - ti(eA x (g e ))(eg e - eA x {g t ))) 

= ^A x (h(eg € ) - h(eA x (g £ )) - H \eA x (g e ))^A x (eg £ - eA x (g e ))) 

= 7fe( A *(M e ^) - h(eA x (g £ ))) = VQ e 

On the other hand 

\ [ l*(e(2wv))(vn x ) + Mdv 

= (l + eA(g e / -i-ii-! (vn x ) + Mdv 

< A(F e ) [ (e 2Mv \-2\w\\v\-l)(vn x )+Mdv 
Jr n 
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since, for each a > 0, the map 

(ae) n - 2 

n>2 



€ ^ 2~ 



= « 2 E 



is increasing. 
Finally 



^ / (w(s,x) ■ v)(v ■ n x ) + F e (s,x,v)dv 
e Jr.* 



a f 

<— I (w(s, x) ■ v ){v ■ n x )+g e (s, x, v)M{v)dv 
e Jr n 
a 

2r> - 



< ^Vg e (s,x) + -C(w)l w(S)X) ^ A(F e )(s,x) 



where 



C(w) := \ I (e^MI^M - 2||iu|| L oo|u| - l){ Vl ) + M(v)dv 
Jr n 



□ 



After integrating in (s,x) both sides of (47), the first term on the 
right hand side of (47) will be absorbed by the Darrozes-Guiraud in- 
formation on the right hand side of (43), so that, with Lemma 4.6, the 
inequality (43) is transformed into 
(49) 

1 „,„ 1 



- 2 H{F e \M w ){t) - —H(F* n \M liew (o : .) t i) 

< — / / V t (s,x)dxds — %r I / VQ t (s,x)dxds 

q Jo Jn 2e J Jqq 



/ // (v — ew(s,x))® 2 :V x w(s,x)F € (s,x,v)dxdvds 
Jo JJnxn N 

(v — ew(s, x)) ■ E(w)(s, x)F e (s, x, v)dvdxds 

C(w) f [ A x (F e )(s,x)dS(x)ds, 
Jo J<9rjnsupp(«)) 



e Jo J Jnxn N 
1 

e Jo J Jnxn N 
ct 



for all t > 0. 

4.5. Control of the outgoing mass flux. In the lemma below, we 
shall prove that the outgoing mass flux A x (F t ) is uniformly bounded 
in Lj oc (H + x <9f2), so that the last term on the right hand side of (49) 
vanishes under the assumption a(e) = o(e) as e — > 0. 
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Lemma 4.7. With the notations of Theorem for each e > 0, each 
renormalized solution F t of the initial-boundary value problem (33)- 
(34)- (35) satisfies the inequality 



(50) 



F e (t,x,v)(v ■ n x ) + dv <-n-^rVQ e (t,x) 

R iV tl{l]) 



+ , 1 / FAv ■ n x fdv 



V / 27r(l - 7]) Jrn 

for all rj E (0, 1), a.e. in (t, x) G R + x dVt. 

Proof. First we recast the Darrozes-Guiraud information in the form 

We = T^UHege) - h{e\ x {g e ))) 
(51) = ^A x (GMG e -G £ - A x (G e ) In A x (G t ) + A x (G e )) 

= * A ^(^ ln (A&)-^ + A ^ ) )' 

since 

A x ((G e - A x (G e )) In A x (G e )) = A x ((G e - A x (G e ))) kvA x (G e ) = . 
Then we consider the integral 

I : = A x G e / (v ■ n x ) 2 + A lMdv 
Jr n 

= A x G t / l|G e /A B G e -l|>Tj(u • n x)l A lMdv 

+ A x G e / l\ Ge /A xGe -i\< v (v ■ n x ) 2 + A lMdv =: h + I 2 
Jr n 

avec 77 g]0, 1[. 

The first term is estimated in terms of the Darrozes-Guiraud infor- 
mation on the boundary, in view of (51): 



R 



iv 



G * 1 



A x G e 



(v ■ n x ) 2 + A lMdv 



- W) L ( G - 1,1 iik) - G - + KG ') iv ■ n ° UMdv 
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As for the second term, one has 

h = A x G e / l\G e /A x G e -i\< v {v • n x )\ A IMdv 
Jn N 

< i / G e l| Ge/A<sGe _i|<„(u • n x f + A IMdv 



1 

1-77 



< / G e (v -n x yMdv . 



>R N 

Putting together both estimates gives 



I = h + h< TT^Vg e (t, x) + / F £ (v n x fdv . 

Km 1 - v Jn N 



Km 1 - rj J r n 

Since 

dz 

is independent of x, we conclude that 



/ (v ■ n x ) 2 + A IMdv = I' z 2 e~ z2/2 
Jr n Jo 



=: J > 



f F £ (t,x,v)(vn x )+dv = -^A x {G e ){t,x) = 

which, together with the previous inequality, leads to the announced 
estimate. □ 



4.6. Convergence to the incompressible Euler equations. At 

this point bring together the relative entropy inequality (43) and the 
boundary control (47), thereby arriving at the estimate 

^H(F e \M w )(t) - ^(FriM^o,),!) 

^ — ZfT / / V £ (s,x)dxds - / / Vg t (s,x)dxds 
e q Jo Jn Ze Jo Jan 

/ II (v — ew(s,x))® 2 : V x w(s,x)F t (s,x,v)dxdvds 

e Jo JJnxn N 

(v — ew(s, x)) ■ E(w)(s, x)F € (s, x, v)dvdxds 



e Jo J In 



t 



+ -C(w) I I A x (F e )(s,x)dS(x)ds. 

e Jo J(9Qnsupp(to) 
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JO Jannsupp(w) Jr" 



Now, statement d) in Theorem 4.1 and the scaling assumption on the 
accommodation parameter, i.e. a(e) = o(e), show that, for each T > 



uniformly in t E [0, T] as e — > 0. 

With (53), the relative entropy inequality (52) is precisely of the 
same form as the inequality stated as Theorem 5 in [29]. One then 
concludes by the same argument as in [29]. 



As recalled above, the convergence of solutions of the boundary value 
problem for the Navier-Stokes equations in the large Reynolds num- 
ber regime is an open problem, as well as the validity of the Prandlt 
equation for the boundary layer. Convergence to solutions of the Eu- 
ler Equation, for general boundary conditions (including the Dirichlet 
boundary condition) are proven under only the most stringent regular- 
ity assumptions. 

On the other hand the onset of von Karman vortex streets and the 
Kolmogorov hypothesis on turbulence based on a non zero energy dis- 
sipation in the large Reynolds (see chapter 5 in [13]) suggests that, in 
general, the limit is not a solution of the Euler equation. This is in 
agreement with Kato's criterion [19] relating the convergence of the 



(53) 




5. Conclusion and final remarks 
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solutions of the Navier-Stokes equations with Dirichlet boundary con- 
dition to a solution of the Euler equations with the vanishing of the 
viscous energy dissipation in a boundary layer with thickness 0(Re _1 ). 

For the Navier slip boundary condition (14) or (13), the inviscid limit 
is established in the present paper is proven under the only assumption 
A ->■ as v ->■ 0. 

The following remarks are in order 

• at variance with previous results [8, 37], no regularity assump- 
tion is required for all the results in the present paper, as only 
estimates are used in the proof; 

• what is proved here is the convergence to a dissipative solution 
(hence to the unique classical whenever it exists); therefore, this 
convergence is also true even if no classical solution exists, or 
even if the L 2 initial data corresponds to a "wild" solution a 
la C. DeLellis and L. Szkelyhidi (see [11]) for which there is no 
uniqueness of the Euler solution; 

• therefore the main goal of the present paper is to show the 
strong similarity between the inviscid limit for the Navier-Stokes 
equations and the fluid dynamic limit for the Boltzmann equa- 
tion; the accommodation coefficient a, the Mach Ma and Strou- 
hal Sh numbers, the Reynolds number Re (see [35], §1.9), and 
the slip coefficient A are related by 

Ma = Sh = e, -L = (l), - = A, 
Re e 

as e — > 0, so that the conditions a = o(e) as e — )• and A — > 
as v — > are consistent. 

In all case the convergence of the solution of the Navier-Stokes equa- 
tions to a classical solution of the Euler equations will imples that 
the energy dissipation vanishes in the limit, as observed by Kato [19]. 
Likewise, the entropic convergence obtained in the present paper in 
the case of renormalized solutions of the Boltzmann equation implies 
that the sum of the entropy dissipation and of the Darrozes-Guiraud 
information at the boundary vanishes with e. 

By analogy with the work of Kato [19], a possible conjecture is that, 
whenever 

a(e) 
lim^ > 0, 

6-s-o e 

the vanishing of both the Darrozes-Guiraud information at the bound- 
ary and of the entropy production implies that the inviscid fluid dy- 
namic limit of the Boltzmann equation is described by a solution of the 
Euler equations. 
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Other, perhaps less delicate problems could be analyzed by the meth- 
ods used in the present paper. For instance, Maxwell's accommodation 
is but one example in a wide class of nonlocal boundary conditions for 
the Boltzmann equation; the Navier-Stokes and Euler hydrodynamic 
limits of the Boltzmann equation should be considered also for such 
boundary conditions (see [35] §1.6). Likewise, the condition (38), which 
may not be verified for all renormalized solutions of the initial bound- 
ary value problem for the Boltzmann equation should be removed on 
principle (as in [29]). 

In other words, there remain many open problems related to the 
issues discussed in the present paper, to which we shall return in future 
publications. 
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